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Abstract 

We study an extended system that without noise shows a monostable dynamics, but when 
submitted to an adequate multiplicative noise, an effective bistable dynamics arise. The stochastic 
resonance between the attractors of the noise-sustained dynamics is investigated theoretically in 
terms of a two-state approximation. The knowledge of the exact nonequilibrium potential allows 
us to obtain the output signal-to-noise ratio. Its maximum is predicted in the symmetric case for 
which both attractors have the same nonequilibrium potential value. 
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I. INTRODUCTION 



During the last few decades a wealth of research results on fluctuations or noise have 
lead us to the recognition that in many situations noise can actually play a constructive 
role that induces new ordering phenomena. Some examples are stochastic resonance in 



zero-dimensional and extended systems 
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IL an; 



a, 



noise-induced transitions jjjj, noise- 
noise- sustained 



induced phase transitions [8J, |9|, noise- induced transport 
patterns 14 . Il5 , Ir]] . noise-induced limit cvcle 17 1, etc. 

The phenomenon of stochastic resonance (SR) — namely, the enhancement of the output 
signal-to-noise ratio (SNR) caused by injection of an optimal amount of noise into a nonlin- 
ear system — stands as a puzzling and promising cooperative effect arising from the interplay 
between deterministic and random dynamics in a nonlinear system. The broad range of 
phenomena — drawn from almost every field in scientific endeavor — for which this mech- 
anism can offer an explanation has been put in evidence by many reviews and conference 
proceedings. See Ref. [l| and references there to scan the state of the art. 

Most of the phenomena that could possibly be explained by SR occur in extended systems: 
for example, diverse experiments were carried out to explore the role of SR in sensory 
and other biological functions [18] or in chemical systems These were, together with 
the possible technological applications, the motivation to many recent studies showing the 
possibility of achieving an enhancement of the system response by means of the coupling of 
several units in what conforms an extended medium P,!^ 



j, lij, L«a, iu, 



l20| , or analyzing the possibility 




of making the system response less dependent on a fine tuning of the noise intensity, as well 
as different ways to control the 

In some previous papers 
in extended systems for the transition 
concept of nonequilibrium potential (23, 



renomenon |21|, |2J 

we have studied the stochastic resonant phenomenon 
aetween two different patterns, and exploiting the 



24j. The nonequilibrium potential is a special Lya- 



punov functional of the associated deterministic system which for nonequilibrium systems 
plays a role similar to that played by a thermodynamic potential in equilibrium thermo- 
dynamics 2^. Such a nonequilibrium potential, closely related to the solution of the time 
independent Fokker-Planck equation of the system, characterizes the global properties of 
the dynamics: that is attractors, relative (or nonlinear) stability of these attractors, height 
of the barriers separating attraction basins, and in addition it allows us to evaluate the 
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transition rates among the different attractors. 

In this work we analyze a new aspect of such a problem studying SR between the attrac- 
tors of the noise-sustained dynamics ^| , that is: the same noise source that induces the 
dynamics, induces the transitions among both structures, and produces the stochastic reso- 
nant phenomenon. Some closely related work correspond to the so called doubly stochastic 
resonance j3, as well as to another previous work j^l related with noise-induced phase 
transitions ja, |9(. In both cases the authors have mainly resorted to a standard mean- field 
approach, or to an estimate of the effective potential, while here we obtain the exact form 
of the noise-induced patterns (stable and unstable ones) as well as the complete form of the 
nonequilibrium potential. In this way we can obtain the transition rates and clearly quantify 
the SR phenomenon by means of the SNR. 

The organization of the paper is as follows. In section|n]we present the model and formal- 
ism to be used. After that, we discuss in section ITTT1 the stochastic resonance phenomenon 
between the homogeneous structure and the inhomogeneous pattern. Finally, we present in 
section HYI some conclusions and future perspectives. 

II. THE MODEL 

We consider a one-dimensional system, limited to the region —L/2 < x < L/2, described 
by the following deterministic equation 

dt(f>(x,t) = d x (D(<f>)d x <f>) + F(<f>), (1) 

assuming Dirichlet boundary conditions (that is 0(±L/2) = 0). This equation can be written 
in a variational form as 



V[<f>] = [ L/2 dx{ - f* dj'DWW) + i[2W)M 2 } ( 3 ) 



where the potential 

I / -,/-■-.. i . . . . / -I- i 

'-L/2 

is a Lyapunov functional (while D(<fi) > 0) for the deterministic dynamics and it is essentially 
the logarithm of the probability density of configuration when Eq.(£Q) is perturbed by an 
additive source of spatiotemporal white noise. 
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The starting point of our stochastic analysis will be Eq. (0) with an additional multi- 
plicative noise, in the Stratonovich interpretation, given by 

where £ is a Gaussian noise with zero mean and correlation (£(x, t)£(x', t')) = 2e5(x—x')5(t — 
t'), being e the noise intensity. For the coefficient of the noise term, g(4>), we adopt 



in order to guarantee that the fluctuation-dissipation relation is fulfilled 

As we are considering the Stratonovich interpretation, the stationary solution of the 
associated Fokker-Planck equation can be written as [28] 

PrfM~exp(-V e ///e), (6) 
where the effective potential V e //[0] is given by 

pL/2 

Veff[<P\=V\<P}-\ dx\nD(<j>). (7) 

J -L/2 

Here A is a renormalized parameter related to e through A = e/(2Ax) in a lattice discretiza- 
tion, where Ax is the lattice parameter 

The extremes of V e ff correspond to the stationary fixed points of the noise-sustained 
dynamics. They can be computed from the first variation of V e ff(4>) respect to <fi equal to 
zero, that is 

nL/2 

5V eff [<t>st\ = - / D(<f>)[d x (D(<j>)d x <j)) + F e// (0)]50(x) dx = 0, (8) 

J-L/2 <t>=<f>st 

where 

F elM = + Aj^W) (9) 

is the effective nonlinearity which drives the dynamics. 

We consider the case of a monostable dynamics in absence of noise 

F = —(f) 3 + b4> 2 , (10) 

and we adopt a model of field-dependent diffusivity which induces an effective bistable 
dynamics. In particular we have chosen 



(D , h > 0), that corresponds to have a larger diffusivity in low density (low 0) regions and 
a lower diffusivity in high density (large 0) ones. With this functional form, F e ff(<j)) in Eq. 
© results 

F eff = -0 3 + 60 2 - ^ = 0(0 - 0x)(0 2 - 0), (12) 

where 0i 2 depend on parameters, in particular on the control parameter A. It is worth 
noting here that in the deterministic problem (A = 0) the reaction term is monostable 
while, as we increase the noise intensity, the effective nonlinear term F e ff becomes bistable 
(within the interval < A < b 2 Do/(8h)) and finally, for A > b 2 Do/(8h) becomes again 
monostable (reentrance effect). Our choice of F and D is one among a plenty of different 
forms for the diffusivity leading to a transition from monostable to bistable and inducing 

28\ . that corresponds exactly to the 



the SR phenomenon (see for instance the one used in 
inverse of the present diffusion coefficient, i.e. -D(0) = D (l + hcf 2 )). Density- dependent 
diffusivities arise in a large variety of systems modelled by reaction-diffusion equations j3| . 
In biology, for instance, population dynamics is usually driven by a diffusivity that depends 
on the local population j30|. We can also find examples in physics, a couple of them are in 
polymer physics (where the diffusion can abruptly drop several orders of magnitude at the 
gelation point js^l) and in diffusion of hydrogen in metals j^. 

A remarkable point is that = is always a root of F e ff = (see Fig. 1). This implies 
(from Eq. (JBJ)) that <p{x) = is an extremum of V e ff[<f>] for all values of A. In what follows 
we will call this structure 0o- 

In order to obtain the non uniform extremes of the potential (and also of the probability 
density) we must (numerically) solve 

±(D(<p st )4-<l>st) +F e// (0 st ) = O, (13) 



dx V dx 

for the stationary regimen profiles <p st (x). This approach allows us to found both, the stable 
and unstable solutions. To analyze their stability we need to calculate 5 2 V e ff, that defines 
a Sturm-Liouville problem, with orthogonality weight D{<p st ). From that analysis it results 
that 0o (defined before) is stable for A > 0, and in the bistability region we have two 
nonhomogeneous symmetric patterns: one unstable <p u (saddle) and one stable S . The 
typical form of these patterns is illustrated in Fig. 2. 

In Fig. 3 we show V e ff[(f) st } vs. A, evaluated on the different stationary solutions. We 
define A c as the value of A at which we have symmetrical stability, i. e. where K//[0o] = 
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FIG. 1: Form of the nonlinearities for the deterministic case (A = 0), bistable case (A = A c ~ 0.8) 
and a monostable case (A = 1.2) in the reentrance region. The vertical scale was changed in the 
deterministic case in order to clarify the figure. The parameters used are: Dq = 1, h = 1/2 and 
6 = 2. Note that = remains as a root in all cases. 



III. STOCHASTIC RESONANCE BETWEEN STRUCTURES 

We are interested in the stochastic resonance phenomena occurring in the above described 
system. For a window of noise intensity the effective dynamics of the system is bistable, 
corresponding to a noise-induced nontrivial dynamics. We will resort to the so-called two- 
state approximation j^, all details about the procedure and the evaluation of the SNR 
could be found in Q]. We consider now that the system is subject, in the adiabatic limit, 
to a time periodic signal of the form b = b + S(t) where S(t) = Ab sm(uot). The usual 
way of rocking the potential is to introduce an additive periodic forcing (or linear periodic 
contribution to the potential). However, in the present case, a small periodic variation of b 
around bo results to be more sensitive to induce the periodic change in the relative stability 
of the two attractors. 

Up to first-order in the amplitude Ab (assumed to be small in order to have a sub- 
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FIG. 2: Fixed points of the noise-sustained dynamics of the problem. We show <f>o = 0, the stable 
homogeneous solution; and both nonhomogeneous patterns: the unstable (saddle) <j) u and the stable 
one 4> s . Here we have A = A c ~ 0.8, while the other parameters values are the same as in Fig. 1 

threshold periodic input) the transition rates Wj take the form 

W\{t) = Hi — cti Ab sm(u t), 

W 2 (t) = // 2 + a 2 A6 fun(u t), (14) 

where the constants //i )2 and ai j2 are obtained from the Kramers-like formula for the tran- 
sition rate 

" det VefM 1 1/2 
J det V eff [(j) u ]\_ 

Here A + is the unstable eigenvalue of the deterministic flux at the relevant saddle point (cf) u ) 
and 



e X p[-(y e// [0 u ]-y e// [^])/e]. (15) 



^1,2 = Wi t2 \s(t)=0 

«i >2 = =F^y !*(*)=*■ (16) 

These results allows us to calculate the autocorrelation function, the power spectrum and 
finally the SNR, that we indicate by R. The details of the calculation were shown in Ref. 
j^j]. For R, and up to the relevant (second) order in the signal amplitude Ab, we obtain 

= 7T (q 2 /Ji +fti/i 2 ) 2 = 7T_ ^1 ^ 2 $ , 17 . 

4a*iA*2 Mi + A*2 4e/ii+/i 2 
7 
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FIG. 3: Nonequilibrium potential V e ff[<p s t], as a function of A, evaluated on the stationary patterns: 
curves correspond to stable (<fi s ), homogeneous (<fo) and unstable (c/> u ) patterns. The arrow indicate 
the point where K//[0o] = V e ff[(j) 8 ], corresponding to A = A c « 0.8. 

where 

$ = / dx D((f)')<f>' 2 d(f)', (18) 

J-L/2 J<t>o 

gives a measure of the spatial coupling strength. In our case 0o = and 

$ = D I"' 2 dx - a rctan[V^(x)] 1 (ig) 

J-L/2 I J 

In Fig. 4 we show the SNR as a function of the parameter A (which is proportional to e). 
The existence of the typical maximum is the characteristic fingerprint of SR. For a window 
of noise intensity values, the system enhances the output to the input periodic signal. We 
see that the maximum SNR occurs at the symmetric situation, that is at A = A c . 

A similar behavior is observed in general for a wide range of values for h and Dq com- 
patible with a bistable effective dynamics. In particular, A c is a monotonically decreasing 
function of h, as we show in Fig. 5. For a given value of h, a numerical analysis of Eq. (JT7jl 
indicates that the maximum of SNR take place at X c (h). Note that, for a given value of A, 
h appears as a additional control parameter that allows a fine tuning of the symmetrical 
condition. Finally, in Fig. 6 we show R c = R(X C ) vs. h in the range of values where Kramer's 
formulae applies 
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FIG. 4: Signal-to-noise ratio vs. A as results from Eq. ([170. Here &o = 2, while other parameters 
remain unchanged. 
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FIG. 5: A c vs. h parameter of diffusivity. For small h values A c , and hence the noise intensity, 
increase monotonically. 

IV. CONCLUSIONS 



The study of SR in extended or coupled systems, motivated by both, some experimental 
results and the technological interest, has recently attracted considerable attention \A y, |4 
fl y, Q- In some previous papers fl fl B 

we have studied the SR phenomenon for the 
transition between two different patterns, exploiting the concept of nonequilibrium potential 
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FIG. 6: Signal-to-noise ratio at A c vs. h. We can see a saturation phenomena as h decreases. 



[23l l24j . In this work we have analyzed the SR phenomenon in an extended system from 
a different point of view, that is studying SR between two attractors of the noise-sustained 
dynamics [3, [3| . 

Some closely related work correspond to the so called doubly stochastic resonance 
as well as to a previous work that is tightly related to noise-induced phase transitions 
[sl, 0] . In both cases the authors have mainly resorted to a standard mean- field approach, or 
to some estimate of the effective potential. Here we adopt a different approach, obtaining 
numerically the exact form of the patterns (both the stable and unstable ones) as well as the 
analytical expression of the nonequilibrium potential. In this way we were able to obtain 
the transition rates and clearly quantify the SR phenomenon by means of the SNR. 

We have seen that the a nonhomogeneous spatial coupling, through density-dependent 
diffusivity, changes the effective dynamics of the system and, in agreement with that 
such nonhomogeneous behavior could contribute to enhance the SR phenomenon. The form 
of the patterns, position of the attractors, barrier's high, explicitly depend on the noise 
intensity. We have found that there are ranges or windows of noise intensities where the 
phenomenon could arise (reentrance) . 

By considering the adiabatic limit and exploiting the two-state approximation we have 
theoretically predicted the occurrence of SR between those patterns. It is worth here re- 
marking that it is the same noise source the one that sustains the bistable dynamics and 
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induces SR for transitions among the corresponding structures. The maximum of the SR 
response occurs in the symmetric case, in agreement with the results found in 0,0]. The 
SR phenomenon is robust respect to variations of the h parameter of diffusivity, and when 
h decreases the SNR maximum increases and shifts toward higher A values. The last fact 
follows from the associated shift of the noise-induced transition to larger noise intensities 
which take place in the spatially uncoupled associated system (i.e. the 0-d system resulting 
from suppressing the gradient term in Eq.(J3J)). 

The consideration of more general forms of couplings in many component systems will 
allow us to analyze SR between noise-induced patterns in activator-inhibitor-like systems. 
We will also study, within the present framework, the competence between local and non- 
local spatial couplings etc. These aspects, together with Monte Carlo simulations of 
the different cases, will be the subject of further work. 
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